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In this paper we show how the method of Zakharov transformations may be used to analyze 
the stationary solutions of the Smoluchowski aggregation equation with a source term for arbitrary 
homogeneous coagulation kernel. The resulting power-law mass distributions are of Kolmogorov 
type in the sense that they carry a constant flux of mass from small masses to large. They are valid 
for masses much larger than the characteristic mass of the source. We derive a "locality criterion", 
expressed in terms of the asymptotic properties of the kernel, that must be satisfied in order for 
the Kolmogorov spectrum to be an admissible solution. Whether a given kernel leads to a gelation 
transition or not can be determined by computing the mass capacity of the Kolmogorov spectrum. 
As an example, we compute the exact stationary state for the family of kernels, K^{mi,m2) = 
(mim2)^'''^ which includes both gelling and non-gelling cases, reproducing the known solution in the 
case ( — 0. Surprisingly, the Kolmogorov constant is the same for all kernels in this family. 

PACS numbers: 82.20.Nk, 82.33.Ln 



I. INTRODUCTION 

Smoluchowski 's coagulation equation provides a mean 
field description of a variety of aggregation phenomena 
P, 2, 3, "4]. The physical picture to bear in mind is one of 
a suspension of particles of varying masses that are mov- 
ing around in d-dimensional space due to some transport 
mechanism. When two particles come into contact they 



stick together with some probability to form a new par- 
ticle whose mass is the sum of the masses of the two con- 
stituent particles. Aggregation is irreversible in the sense 
that large aggregates are not permitted to break up into 
smaller ones. If it is assumed that there are no spatial 
correlations between aggregates then the concentration 
of particles of mass to, c(m, t), obeys the Smoluchowski 
kinetic equation: 



dc{m, t) A 



dt 
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The kernel K{mi, TO2, to) and the constant A control the 
rate at which particles of masses toi and TO2 react to 
create particles of mass to = toi -|- TO2. The ^(to — TOq) 
term provides a source of particles of mass TOq such that 
the total rate of mass input is given by Jq which we take 
to be constant in time. J[c] represents the mass flux 
which is functionally dependent on the entire spectrum, 



c(rn,t). Thus the S{m — M) term provides a sink by 
removing particles from the system whose masses exceed 
M. 

The kernel must be symmetric in its first two argu- 
ments, K{mi,m2, m) = K{m2, toi, to), if it is to describe 
a physical aggregation process. Owing to the presence of 
the delta functions, the kernel is effectively a function of 
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two arguments rather than three and is usually written as 
-R'(mi, TO2). We include the explicit dependence on the 
third argument only for notational convenience. After 
writing K{mi,m2,m) = K(mi,'m2), some simple ma- 
nipulations reduce Eq. |^ to the more "standard" form 
often considered in the literature: 



dc{m,t) A 



dmidm2K{mi, TO2)c(mi, t)c{m2, t) 



[5{m — mi — 777,2) — ^("7- — mi) — 5{m — m2)] 

+ —i{m-mQ)-'^3(m-M). (2) 
mo M 

Note that the addition of the aforementioned source and 
sink terms allows a time independent steady state to be 
reached in the limit of large time. This is the main sub- 
ject of this paper. We shall be concerned with the situa- 
tion where mo and M —>■ 00, bearing in mind that 
the presence of a sink at infinity may be required even at 
finite times in the case of the so-called "gelling" kernels. 
As mentioned already, we only consider here sources for 
which the total flux of mass into the system, Jo, is con- 
stant. In the turbulence literature the interval of masses 
for which mo m <^ M is called an inertial range. The 
stationary states considered in this article are valid in 
this range. 

The details of the transport mechanism and stick- 
ing probability are assumed to be built into the kernel, 
K(mi, m2), of the Smoluchowski equation. Different ker- 
nels arise in different physical contexts and determine 
how the solution of the equation should be interpreted 
physically. We refer to Refs. 0, 0| for a short list of 
commonly considered kernels and their physical and/or 
mathematical contexts. 

Most of the kernels of physical interest are homoge- 
neous functions of their arguments. We shall denote the 
degree of homogeneity of the kernel by (. That is 



K(hmi, hm2, hm) — h'^Kijni, m2, m). 



(3) 



This homogeneity need not be uniformly weighted be- 
tween the two arguments. Following P|, we introduce 
exponents, ^ and u to take into account this fact: 



A'(mi, m2, m) ^ m^m2 form2 3>mi. 



(4) 



The exponents /i and v satisfy ^ + 11 — C,. Let us consider 
a couple of simple examples to clarify our notation. The 



kernel K{mi,m2) = X{m^ 



l+e 



has C 
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fi — and = 1 + e whereas the kernel K [mi, 1112) 



1/3 



1/3 



)(mi 



-1/3 



m. 



-1/3 



) has A = 0, /X 



-1/3 



and 7/ = 1/3. These basic properties of the kernel are all 
we shall require for what follows. 

In this paper, we study the steady state behavior of 
c(m) when mo ^ m <C M. In Sec.m using dimensional 
analysis, we derive the large mass dependence of c{m). It 
is also shown that the power law spectrum corresponds 
to a constant flux of mass in mass space. In Sec. IIIII we 
show that the Smoluchowski equation is mathematically 



very similar to the kinetic equation for 3- wave turbulence. 
Using Zakharov transformations from 3- wave turbulence, 
we rederive the mass spectrum as well as compute the 
amplitude also known as the Kolmogorov constant. The 
characteristic mass of the source, mo, does not appear in 
the dimensional argument. In Sec. IIVI we find the con- 
ditions under which this assertion is correct when we ad- 
dress the question of the locality of the mass cascade. In 
Sec.0 we discuss the notion of mass capacity of the Kol- 
mogorov spectrum and show how it distinguishes between 
gelling and non-gcUing kernels. In Sec. IVII we explicitly 
compute the Kolmogorov spectrum for a one-parameter 
family of kernels given by K,^{mi, 1122) = (mim2)^/^. We 
find that the value of the Kolmogorov constant is the 
same for all models in this family. Finally, we end with 
a summary in Sec. IVIII 



II. DIMENSIONAL DERIVATION OF THE 
STATIONARY SPECTRUM 

Before proceeding into detailed analysis of the station- 
ary states of model Eq. |^ let us first describe intuitively 
what we mean by a Kolmogorov solution by employing 
a simple scaling argument. We shall use the simplified 
form Eq. jSJl for brevity. The stationary energy distribu- 
tion of forced hydrodynamic turbulence is described by 
the famous Kolmogorov 5/3 spectrum (for instance, see 
5]). This spectrum, postulated from dimensional con- 
siderations, carries a constant flux of energy from large 
scales to small by means of vortex- vortex interactions. 
The analogous cascade for the Smoluchowski equation is 
a cascade of mass from small particles to large mediated 
by the coagulation of aggregates. The Kolmogorov spec- 
trum for aggregation carries a constant flux of mass. 

The physical dimensions of the various quantities ap- 
pearing in Eq. ^ are as follows : [c] = M^^ L ''i 
[J] = ML-'^T-i and [A] = M-'^L'^T-^. If we now take 
the combination c ~ J''X"m^ , simple dimensional anal- 
ysis requires that we choose 7 = 1/2, a — —1/2 and 
f3 = — (C -I- 3)/2. Dimensional considerations therefore 
lead us to a spectrum of the form 



elm 



(5) 



The characteristic mass of the source, mo, does not ap- 
pear in our dimensional argument on the basis that we 
expect this solution to be valid for masses much greater 
than mo. We find the conditions under which this asser- 
tion is correct in Sec. IIVI when we address the question 
of the locality of the mass cascade. 

The exponent Eq. |(SJ) is not new. It appeared in early 
work by Hendriks, Ernst and Ziff as the scaling of 
the post-gel stage of gelling systems. Their work makes 
an implicit connection between this scaling and the fact 
that there is a mass fiux out of the system in the post-gel 
stage. It was then derived explicitly for the Smoluchowski 
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equation with source term by Hayakawa for a particu- 
lar family of kernels but without making any connection 
with the physical role played by the mass flux. 

That the spectrum Eq. jSJ corresponds to a constant 
flux of mass in mass space is easily seen from the following 
scaling argument. We express the local conservation of 
mass by means of the continuity equation 



dmc{m, t) dJ{m, t) 



dt 



dm 



(6) 



where 

dJ{'m,t) mX 



dm 



= / dmidm2[K{mi,m2)c{mi,t) 

2 Jo 



c(m2, t) {S{m — mi — m2) — S{m — mi) — d{m — m2)} |7) 

We now assume a stationary spectrum, c{m) = Cm~^ . 
By introducing new variables, mi — m/ii, m2 = m^2 
and using the scaling properties of the kernel we deduce 
that 



dJ 
dr 



oc m 



2+C-'2x 



(8) 



with the constant of proportionality being given by the 
integral expression which remains after scaling out the m 
dependence of the RHS of Eq. . Thus 



J(m) cx 7713+^-2^. 



(9) 



It is clear from Eq. ^ that in order to have a stationary 
state, J must be independent of m which determines the 



exponent of the Kolmogorov spectrum as 



XK 



3 + C 



(10) 



Of course we cannot determine the Kolmogorov constant, 
C, from such a scaling argument. In addition the validity 
of our scaling argument depends on the convergence of 
the various integral expressions which have been hidden 
behind proportionality signs. 

In this paper we address these short-comings by com- 
puting the exact stationary solutions of Eq. ||2J) using 
the method of Zakharov transformations borrowed from 
the theory of wave turbulence. We obtain the exponent, 
xk expected from scaling considerations and the value 
of the Kolmogorov constant, C. An answer is obtained 
for arbitrary homogeneous kernels. However the analysis 
involves the exchange of orders of integration on the RHS 
of Eq. |5J . It is thus necessary to check a posteriori that 
the RHS is convergent on the prospective spectrum in 
order that it be an admissible solution. This check leads 
to a "locality criterion" , namely 



^ - + 1 > 0, 



(11) 



which must be satisfied by the kernel in order that the 
Kolmogorov spectrum be realizable. 

III. ZAKHAROV TRANSFORMATION FOR 
SMOLUCHOWSKI EQUATION 

To find the stationary solutions of Eq. 1^ in the situ- 
ation TO ^ Too, Jo = constant, we must solve: 



A 



— / dmidm2 K{mi, TO2, to) c{mi)c{m2) (5(m — TOi — TO2) / dmidm2 K{mi, to, TO2) c(to)c(toi) J(m2 — to — TOi) 



dmidm2 K{m, TO2, toi) c(TO)c(m2) S{mi — TO2 — to). 



(12) 



r 



Structurally this equation is very similar to the kinetic 
equation for wave turbulence with a 3-wave interaction. 
For an introduction to the theory of wave turbulence see 
. A useful trick for finding the stationary power law so- 
lutions of such equations was devised by Zakharov ^ 
in the late 60's and is easily applied here. Restricting our- 
selves to power law solutions of the form c(to) = Cto~^, 
we apply the following changes of variables, 



(toi,TO2) 



111^ III 
~i ' T ) 

TO, Too 



TO TOTO. 
(TOi,TO2) -> ( — , - 

m\ m\ 



(13) 
(14) 



to the second and third integrals in Eq. respectively. 
Dropping the primes on the transformed variables and 



using the homogeneity and symmetry properties of the 
kernel we obtain the equation 



= 



XC' 



dmidm2[K{mi,m2,m){mim2) ^m 



2rE-C-2 



(5(to — TOl — TO2X]l5) 



It is immediately evident that the integrand is identically 
zero for 2x — — 2 = 1 from which we get the same 
Kolmogorov exponent. 



c 



(16) 



obtained in Sec.^by a scaling argument. The value of the 
Kolmogorov constant can be determined by considering 
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the local mass flux defined from Eqs. © and Re- 
stricting ourselves to spectra of the form c(m) = Cm^^, 



the Zakharov transformation allows us to write Eq. Q 
in the form 



dJ{m,t) _ XmC^ 



dm 



dmidm2K(mi,m2,m)(mim2) 



2X-C-2 2a:-C-2\ \ 

I — m2 I o(m — 77ii — 7712) 



r 



where 

lix) 



i-^r'-2-Mr^X'5(i-Mi-A*2)](i7) 



From this we deduce that the flux is given by 

= ^^l!l(^^3+C-2x, 
^ ' 3 + C-'2x 



(18) 



In the steady state, x = xk = (3 + 0/2 and the flux 
must be a constant equal to Jq. Thus we have 



\C'^I{x) _s,+c~2x 



x-^xk 3 + — 2a; 



(19) 



We know that I{xk) = so we must apply I'Hopital's 
rule to evaluate the limit to arrive at 



, _ A 2 dl 



(20) 



and hence 



C 



'2 Jo dl 
X dx 



(21) 



The Kolmogorov solution is therefore, 

c(m) — Cm~^^ , 
with C given by Eq. ^ and Xk given by Eq. 



(22) 



The support of the integrand on the RHS of Eq. ^ 
is shown in Fig. ^ Since the kernel and the mass dis- 
tributions which we study are scale invariant the only 
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FIG. 1: Support of the integrand of Eq. Q 



possible sources of divergences are at infinity and at the 
two points (0,m) and (m,0) where the contour of inte- 
gration intersects the axes. Let us now study carefully 
the behavior of the integrand near these points for power 
law mass distributions. 

The behavior at infinity is easy. As m2 — s- cxd along the 
lower contour the integrand looks like 



IV. LOCALITY OF THE MASS CASCADE 

In the analysis of the previous section we have freely 
split the integrand on the RHS of Eq. |^ and exchanged 
orders of integration to derive the Kolmogorov spectrum. 
In order to justify these manipulations we must demon- 
strate a posteriori that the original collision integral is 
convergent on the Kolmogorov spectrum. We do this in 
this section. 



/^(m, m2, m + m2)c(m)c(m2) ~ ^- (^3) 



The integral is therefore convergent as rn2 — > 00 if 

v-x<-\^x>v^\. (24) 

The same criterion is obtained along the upper contour. 
The convergence near zero requires a little care. Near 
TOi = the integrand looks like 
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c(toi) [K{mi, m — mi, m)c{m — mi) — K{mi, m, m + TOi)c(m)] 
c{mi) 



d d 
™i ( 7^ {K{mi,£„m)c(i)\^^- c{m)— {K{mi,m,C))\^^^ 



+ o{m\) 



m 



1— 



Note the cancellation of the leading order terms in the 
Taylor expansion on the second line above. The corre- 
sponding integral is convergent as toi — > for 

^ + l-a;>-l^.i:<^ + 2. (25) 

The same criterion is obtained if we look near m2 = 0. 
We conclude from Eqs. 124|l and (|25|l that a power law 
mass distribution, Cm^^, yields a convergent collision 
integral if x lies in the interval [v + 1, /i + 2]. The exis- 
tence of such an interval of convergence puts a constraint 
on the asymptotic behavior of the kernel, namely 

^-i^+l>0. (26) 

We now must address the question of when the Kol- 
mogorov spectrum derived in Sec. IIIII lies in this inter- 
val of convergence. The answer is surprisingly simple. 
Remembering that /i + = C, it is immediately evident 
from Eq. (|16|l that the Kolmogorov spectrum lies midway 
between the two constraints Eqs. H24(l and H25|l . There- 
fore if an interval of convergence exists for a given kernel, 
the corresponding Kolmogorov spectrum is an admissible 
stationary solution of Eq. and it lies at the midpoint 
of the interval of convergence. 

We call Eq. I|2H(I a locality criterion since systems for 
which it is satisfied can be characterized in the station- 
ary state by a local mass flux, Jq. When the spectrum 
is local, the details of how we take the limits mo — > 
and M ^ cx) to produce a large inertial range are in- 
consequential since all integrals converge. If Eq. H26|l is 
not satisfied then presumably the final stationary state 
depends on the details of the source/sink and is therefore 
non-universal. 

We note that the kernel m^^'^-l-mj^'^, mentioned in the 
introduction, is marginal in the sense that it violates the 
locality criterion for any finite e. It has been shown [llj 
that this kernel undergoes instantaneous gelation so per- 
haps there is some connection between this phenomenon 
and the concept of locality. In addition, the generalised 
sum kernel, K{mi, 1112) = ^ + m,2^ , which violates the 
locality condition for u> 1, was studied extensively by 
Krapivsky et al. p^ . Il3| |. They found that in this case, 
the system does not reach a steady state but rather con- 
tinues to evolve very slowly on a logarithmic timescale 
for all time. 

In closing this section it should be noted that a rig- 
orous understanding of the conditions under which the 
stationary state depends only on the local flux is one of 
the missing pieces in the theory of hydrodynamic turbu- 
lence. 



V. FINITE AND INFINITE CAPACITY CASES 
- GELLING AND NON-GELLING KERNELS 

It was found in the 60 's [3 that the solution of Eq. (O 
for certain kernels violates mass conservation within a fi- 
nite time, t*. When this violation occurs, lim,„^oo P{m) 
becomes finite. In the late 70's it was found that mean- 
ingful solutions exist post-i* and the violation of mass 
conservation was given a physical interpretation in terms 
of what is now termed a "gelation transition" 0,0,^3- 
Gelation occurs when there is a finite flow of mass to 
an infinite mass cluster ("gel"). As a consequence, the 
total mass of the normal ("sol") particles is no longer 
conserved. In order to avoid inconsistencies, the gel par- 
ticles must be considered as those clusters whose mass 
diverges as the size of a finite system is taken to infin- 
ity. It is now well known 0] that the solutions of Eq. (0) 
undergo gelation for kernels having C > 1- 

The gelation criterion, C > 1; can be given a very sim- 
ple physical interpretation by examining the mass capac- 
ity of the Kolmogorov spectrum. If we continue to add 
mass to the system at a constant rate Jq then we know 
that the final steady state is given by the Kolmogorov 
spectrum Eq. H22|l . When the total mass contained in 
this solution is finite then mass conservation must be vi- 
olated at some time since the total mass supplied to the 
system grows linearly in time. The total mass capacity 
of the Kolmogorov spectrum is finite when 




This integral in convergent at its upper limit when 1 — (3+ 
0/2 < — 1 or ^ > 1. Thus gelation can be seen as a kind 
of safety valve which allows mass to fiow out of the system 
when the Kolmogorov spectrum is incapable of absorbing 
all of the mass supplied to the system. Conversely, one 
would expect intuitively that infinite capacity systems 
should not exhibit gelation. 

VI. EXAMPLE : THE FAMILY OF KERNELS, 

In this section we explicitly evaluate the Kolmogorov 
constant, C, for the family of kernels, if^(mi,m2,m) = 
(mim2)'»/^. These kernels have ^ = v = C/2 so that 
the corresponding Kolmogorov spectrum always satisfies 
the locality criterion, Eq. (|26|l . The family includes both 
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gelling and non-gelling kernels. In general to compute, C, 
we need to evaluate the following integral at a; = (3 + C)/2 



dl 

dx 



dfiiKifii, 1 - ^1, 1) - fi,))-- [-2^^ logMi - 2(1 - log(l - fi,) 



■ iog(/ii(i - Mi))(i - - (1 - 

r 



(27) 



This is obtained from Eq. p7|l by integrating out ^2 
and differentiating with respect to x. When we set 
K{jjLi,^2, 1) — (a^iA^2)''^^ and x — {?> + C)/2 in this ex- 
pression we find, rather surprisingly, that all dependence 
on C cancels out and we are left with 



dl_ 

dx 



(3+f)/2 



\ MilogMi + (l-Mi)log(l-/^i) 
"Ml 



= 47r 



3/2/1 

Ml (1 

(Mathematica) 



(28) 



Hence the Kolmogorov solution for all kernels in this fam- 
ily is 



^ ' 27rA 



(29) 



To close we note that we can check our answer indepen- 
dently for at least one case. For the constant kernel with 
zero initial concentration, an exact solution of Eq. (0) 
has been known for some time. The details can be found 
in 18] . This solution is 



Cm.{t) = ^ Ck{t) S{m - kmo) 



(30) 



k=l 



with 



Ck{t) 



]=co 



(2j + ifmoTT^ 



2AJoi2 



(31) 

The < ^ 00 limit of this expression exists can be calcu- 
lated by replacing the sum by an integral in the limit of 
large t. This integral can be expressed in terms of gamma 
functions. One finds 



lim Cfe(i) 

t — >oo 



Jo r(fc - 1) 

27rAmo r(fc + 1) 

-k^i for k :$> I. 



Jo 



2nXmo 

Setting TTiQ = 1 we recover the result of our earlier com- 
putation of the Kolmogorov spectrum for m 3> 1. 

As pointed out to us by one of our referees, the con- 
stant, C, has also been computed 0,0] for the gener- 
alised sum kernel 



K{mi,m2) 



(32) 



We computed the integral H27|) for this kernel using Math- 
ematica and found the Kolmogorov constant to be 



C 



^Jo(l-e)cos^ 
4A7r 



(33) 



as found in using completely different methods. 



VII. CONCLUSION 

To summarize, we have shown how the notion of a 
mass cascade analogous to the Kolmogorov energy cas- 
cade of hydrodynamic turbulence is relevant to under- 
standing the stationary state of the Smoluchowski equa- 
tion with constant mass production term. Furthermore 
we have shown how the exact stationary spectrum may 
be computed using the method of Zakharov transforma- 
tions and given some criteria for assessing the physical 
validity of this solution. We have not made any at- 
tempt to address the important question of the valid- 
ity of the Smoluchowski equation itself in describing the 
statistics of particular aggregation models. The mean 
field assumption leading to Eq. |^ can be violated in 
several ways as discussed in [l^ . Of particular relevance 
to lattice aggregation models is the case where fluctu- 
ations dominate the statistics and invalidate the mean 
field Smoluchowski equation [T^ I2(tI | . In a future publi- 
cation |2l| we shall address this issue for the particular 
case of constant kernel stochastic aggregation where the 
presence of fluctuations leads to a renormalization of the 
constant A. The techniques developed in this paper will 
allow us to find the renormalized Kolmogorov spectrum 
as the stationary solution of a modified Smoluchowski 
equation. 
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Note added in proof lems [22 for a class of kernels arising from the study of 

galactic mergers in astrophysics. 

We have found that V.M. Kontorovich has recently ap- 
plied the Zakharov Transformation to aggregation prob- 
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